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Q Abstract 

We study the dynamics of non interacting thermal atoms embedded in structured optical lattices 
^ with non trivial geometry. The lattice would be generated by two counter propagating modes with 

parabolic cylindrical symmetry and we concentrate on the quasi conservative red detuned far-ofF- 
^ resonance regime. The system exhibits quasi periodic and chaotic behaviors whose probability can 

a 

^ be controlled by varying the intensity of the beams. The spectral density of the trajectories is used 

as a chaos signature. An analysis of permanency times for chaotic trajectories that visit more than 
J> one potential well reveals a distribution with a long tail. 
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I. INTRODUCTION 



The advent of optical cooling and trapping of atoms gave rise to the experimental and 
theoretical search of classical and quantum chaotic effects in the evolution of atomic center 
of mass motion [THH]. It was found that the spatial and temporal stochastic character of 
spontaneous emission can induce random walks on the atoms during laser cooling processes 
[1]. Anomalous transport properties were discovered, like the fact that both variance and 
mean time for atoms to leave the laser beams could become infinite; the atoms movement 
became then a particular class of Levy flights [5]. 

Other interesting effects arise in optical lattices built from temporally modulated standing 
waves made using two Gaussian counter propagating laser beams. If the modulation is chosen 
to generate an effective periodically driven rotor, the classical dynamics can become chaotic 
and a quantum manifestation is the localization in momentum space of the ultra cold atoms 
[2] . In Ref . [3] the quantum dynamical tunneling of ultra cold atoms between classical islands 
of stability was also reported; this shows that the presence of chaos in phase-space can yield 
substantial consequences on the quantum tunneling rate between classical regular regions, 
an effect that has since become known as chaos-assisted tunneling [6]. 

Spatial light modulation to study chaotic dynamics of cold atoms has also been imple- 
mented in the context of optical billiards. To create an optical billiard, a laser beam can 
be deflected at different angles synchronously so that an arbitrary two-dimensional light 
pattern is formed in a plane perpendicular to the optical axis. An additional standing wave 
aligned perpendicular to the billiard plane confines the atomic motion to two dimensions 
[?]. This scheme has been used to investigate the chaotic and regular dynamics of atoms 
on well-known billiards fS^j "9]. A technique named echo spectroscopy has been developed to 
study quantum coherence in the evolving atomic system [10]. Introducing variable Gaussian 
beams waists into this scheme simulates soft- wall billiards 

In this work we study the semiclassical dynamics of non interacting cold atoms in optical 
lattices built from two counter propagating beams with transverse structure. We numerically 
show that, even in the simplest case of far-off- resonance quasi conservative dynamics, the 
transversal geometry of the beam can be used to generate classical chaotic dynamics. This 
occurs both in the tight binding regime (similar to an optical billiard) and in the case 
that initial conditions allow the atomic transport between different potential wells. All the 
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reported results will consider optical lattices with parabolic cylindrical geometry although 
other geometries could also yield similar results. 

In general, the dynamics of cold atoms in optical lattices is highly dependent on the optical 
and atomic parameters that nevertheless are feasible to be controlled. Dilute atomic samples 
with predetermined cold atom-cold atom interactions can be used to probe experimentally 
single-atom and many-atoms phenomena. The detuning between the light frequency and 
the chosen atomic transition frequency is used to regulate both the relevance of dissipative 
effects, related to spontaneous emission, and the depth and sign of the effective light potential 
affecting the atoms. 

It is important to emphasize that structured electromagnetic beams are experimentally 
feasible [12], with potential applications for the manipulation of cold atomic systems in 
the semiclassical and quantum regimes[T3]; some of these applications have already been 
implemented [HI [T3] . 

II. OPTICAL LATTICES WITH TRANSVERSE PARABOLIC SYMMETRY 

The electromagnetic (EM) waves with cylindrical symmetry and transverse structure 
exhibit many interesting features. Ideally, these waves have an intensity pattern invariant 
under propagation along a given axis that we take as the z axis. The best known examples 
correspond to Hermite and Bessel waves. The first have Cartesian symmetry and the latter 
have circular transverse symmetry [JL6J. Elliptic symmetrical waves are known as Mathieu 
waves [TTj. The fourth and last separable example, has a transverse structure naturally 
described in terms of parabolic cylindrical coordinates (m, f , z) [18] 

X \ % y = -{u \ % vf ^ z = z (1) 

where x, y and z are the Cartesian coordinates, and u G {—oo, oo) and v G [0, oo). Surfaces 
of constant u form half confocal parabolic cylinders that open towards the negative x axis, 
while the surfaces of constant v form confocal parabolic cylinders that open in the opposite 
direction. The foci of all these parabolic cylinders are located at x=0 and y=0 for each z 
value. The separability of the wave equation in such a coordinate system, allows writing the 
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electric and magnetic fields of a monochromatic mode as [19] 



= -^f ^)9,,Mvl/, - A^™^d,,m., (2) 

J —n 

^ia In I tan 'p/2\ 

2le(a;(/?) 



2to(a;v2) 



2y/7r sin ip 



A^r^ and ^1™) are proportional to the amplitude of the transverse electric and transverse 
magnetic modes, and \E'k is a solution of the scalar wave equation with k representing the 
whole set of numbers that specify a particular mode. That is, the parity p, the wave vector 
component along the main axis of propagation kz, the frequency uj = C\/k\ + k\ and the 
separation constant a involved in a specific solution of the scalar wave equation. 

In this work we shall report results for parabolic cylinder electromagnetic modes of order 
a = 0. Their intensity and polarization structure are illustrated in Fig. [TJ These modes, 
for both even and odd parities p, have been experimentally generated by means of a thin 
annular slit modulated by the above angular spectra 2lp(a; y?) [SU] . 

An optical lattice with parabolic cylinder structure can be generated by the superposition 
of two counter propagating beams with that symmetry. For simplicity we restrict our study 
to TE modes with even parity. The structure of the lattice along the cylindrical symmetry 
z-axis will be that of a sinusoidal standing wave. 



III. SEMICLASSICAL FORCE ON THERMAL NEUTRAL ATOMS 

Under standard conditions, the interaction between a two level atom and an EM wave 
with a frequency close to resonance has an electric dipole nature with a coupling factor 

9^ = (^i2i^x ± idy) ■ E. 
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Figure 1: (a) Transverse intensity and (b) electric field of a TE electromagnetic wave with 

parabolic cylindrical symmetry. The wave parameters are kz = 0.975u/c, a = and the 
associated scalar field has even parity. Length is measured in units of the light wavelength 



In a semiclassical treatment, the gradient of g, 

Vg = {a + i(3)g, 



defines the force experienced by the atom. The expression for the average velocity dependent 
force, valid for both propagating and standing beams, is 



(f ) = hfp' [(v • a 



1 — p 

'l+p 



+ r/2]/3+ [(v/3) 



In this expression 



f = r/[r(i + p') + 2v ■ a[i - p/p' - p] [p'/{i + p)]], 



(3) 



(4) 



r = 4fc^|/xi2p/3/i is the Einstein coefficient, 6uj = u — uq denotes the detuning between the 
wave frequency u and the transition frequency ojq, p = 2|(7p((r/2)^ + Su"^) is known as the 
saturation parameter, linked to the difference D between the populations of the two levels 



of the atom, D = + p), and finally p' = 2\g\ 
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with 7' = (v ■ a.){l -p){l +p)~^ + 



T/2 + i[{^r ■ f3) -6cu]. 

The dissipative term (v ■ associated with a Doppler shift, as well as other velocity 
dependent terms in Eq. ^ are expected to be small for slow atoms, particularly within the 
red detuned far-off- resonance regime [23] . In such a case Eq. ^ takes the simpler expression 



(f) = h: 



p 



p 

5 



-f3 — 6ucx 



(5) 



However, we keep the velocity dependent terms in our numerical calculations in order to 
prevent disregarding potentially relevant effects. 

We consider that the cylindrical symmetry axis z of the light field configuration is oriented 
along the vertical direction and gravity force is included. 

IV. ATOM SEMICLASSICAL TRAJECTORIES 

The numerical simulations [24J consider a TE laser beam detuned 67 nm to the red of 
the 5 '^S\i2 - 5 transition at 795 nm of ^^Rb with irradiance in the range ~ 0.5 — 22 

kW/ cm^. As natural length and time units, we take the laser wavelength and the inverse 
of the Einstein coefficient F, which is 3.7 x 10''s~^ for the state 5 ^Pi/2 of '^^Rb. 

Given an irradiance, initial conditions for a cloud of a hundred non interacting atoms are 
generated randomly within a circle of radii 20A centered on the axis of the beam. The random 
velocities correspond to a temperature in the range of 2.9-3.1 /iK for the movement in the 
plane XF, perpendicular to the beam, and about 0.2 /iK along the z axis. The anisotropy 
in the initial velocities enhances the effects of the nontrivial transverse structure since it 
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Figure 2: Percentage of non trapping (circles), quasiperiodic (stars) and chaotic trapping 
(squares) trajectories as a function of the light beam irradiance. The initial conditions for 
the atomic cloud correspond to random initial positions within a circle of radii 20A 
centered on the axis of the beam and velocities corresponding to a temperature in the 
range of 2.9-3.1 /iK for the movement in the plane XF, perpendicular to the beam, and 

about 0.2 /iK along the z axis. 
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Figure 3: Three illustrative examples of quasiperiodic trajectories of a cold atom trapped 
in a parabolic cylinder lattice: (a) the movement takes place within the second lobe to the 

left of the symmetry axis (see Fig. [T]); (b) the movement is confined to a high intensity 
lobe closest to the symmetry axis; (c) a typical quasiperiodic trajectory that involves more 
than one optical lobe. For clarity the time interval plotted is a few percentage (~ 10%) of 
the 1.25 X 10*^ time considered in the numerical calculation. Length is measured in 

units of the light wavelength A. 

"freezes" the z-direction degree of freedom. A mechanism for achieving initial anisotropic 
velocities at low temperatures could be based in the previous use of a dissipative optical 
lattice with an anisotropic modulation of the laser-atom interaction parameters [25|. 

In Figure [2] we illustrate the proportion of non trapped, regular trapped and chaotic 
trapped trajectories as a function of the irradiance. For instance for an irradiance of ~4.5 
kW/cm^ about 27% of the atoms are not trapped by the lattice and escape with their route 
either almost confined in a transverse plane to the beam or with their velocity almost parallel 
to it. The remaining atoms exhibit a trapped chaotic (28%) or quasi periodic motion (45%). 
An atom will be considered trapped if for times lower than 1.25 x lO^F"^ it remains at a 
distance to the axis lower than 80A while, in the axial direction, it remains within at most 
5A of its initial z value. An irradiance of ~ 2 kW/cm^ is required for trapping ~ 50% of 
the atoms. Notice however, that for an irradiance as low as ~0.5 kW/cm^ there is a high 
probability of observing chaotic trajectories within the trapped ones. For an irradiance of 
~ 6 kW/cm^ ~ 82% of the atoms are trapped by the lattice, about ~ 35% have a chaotic 
motion and ~ 47% a regular one. In the following, the cases that illustrate the atoms 
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dynamics refer to the latter irradiance. 

In most of the regular trajectories, the atoms are confined to a single lobe of the lattice 
(Fig. [3^, ^p) although less than 10% of the observed regular trajectories can involve two 
lobes (Fig. ^) or more than one well along the 2;-direction. We could not find regular 
trajectories visiting more than two lobes in the XY plane. The quasiperiodic character of 
the motion can be easily verified by evaluating the power spectra of any of the components 
of the position or velocity vectors as a function of time as shown in Fig. |4j 

As for chaotic trajectories, for the initial conditions we considered, they give rise to 
confined motion within one wavelength in the z direction and transversal motion in either 
one, two or several lobes in the labyrinths created by the optical lattices in the XY plane, 
as illustrated in Fig. [5j In the most common chaotic trajectories, the atom remains in a lobe 
for certain time and then goes to another lobe of the lattice in a very irregular way that 
practically covers the accessible configuration space in the optical lattice. Fig. |5]3. 

For atom trajectories involving many lobes, the lobes are visited in an irregular order 
and an atom can visit the same lobe several times. For Fig. the atom visits most of the 
lobes in any interval (to)^o + ^t) of duration At = 1.25 x lO'^F"^ irrespective of the initial 
time to- 

In Fig. |6] the power spectra of the x coordinate of the trajectories shown in Fig. [5] is 
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Figure 4: Power spectra of the coordinate x{t) for the three quasiperiodic trajectories 
shown in Fig. [3] showing the typical well defined isolated periodicities expected for a 
regular dynamics. Only the low frequency end of the spectra is shown. The frequency is 

measured cycles per unit time F"^. 
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(a) (b) (c) 

Figure 5: Three illustrative examples of chaotic trajectories where: (a) the movement takes 

place within the second lobe to the right of the symmetry axis (see Fig.g; (b) the 
movement is confined to the two highest intensity lobes closest to the symmetry axis; (c) 
the trajectory practically covers the accessible configuration space in the optical lattice. 
For clarity the time interval plotted is a few percentage (~ 10%) of the 1.25 x 10^ F"^ 
time considered in the numerical calculation. Length is measured in units of the light 

wavelength A. 



presented. The broad band structure of those spectra is a signature of the trajectories 
chaotic character. 

Looking for other signatures of chaos, the phase space trajectory can also be also analyzed. 
In Fig. [7]we illustrate such trajectories in {y,Vy) space. Both Fig. [sj^ and Fig. [T] point to a 
dense covering of the accessible phase-space established by the light labyrinth and the atom 
initial conditions. 

In Fig. |8| we illustrate the distribution of the permanency time Tperm within a lobe in 
a log-log plot for the trajectory showed in Fig. [Sjo. The sampling was taken considering a 
total evolution time Tp = 2.5 x lO^F^^ in which the atom transits 82951 times from one 
to the other lobe. Analysis with shorter total evolution time (we studied Tp in the interval 
[0.5 X 10^r^^,2.5 X lO^F"^]) gave frequency distributions of Tperm with the same structure; 
that is (i)the same characteristic minimum time an atom expends within a lobe (a time 
naturally determined by the initial kinetic energy); (ii) local very stepped maxima followed 
by strong decays which are then followed by the next maxima; (iii) the longest permanency 
time Tiongest has a very small frequency; (iv) if the longest permanency time is not considered 
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the distribution of the local maxima is approximately linear, which would then lead to a 
power law statistics for that maxima; (v) Tiongest increases as Tp increases. This indicates 
to us that the momenta of the distribution is not well defined, like in a Levy distribution. 

Partial trajectories with nearby values of the permanency time Tperm exhibit similarities 
as illustrated in Figs. [9] and 10 Fig. refers to the trajectories with Tperm equal to 



the minimum time an atom must expend in a lobe. For them, the atom enters through 
one of the potential maxima and leaves through the other one. In most cases the partial 
trajectories associated with the local maxima of the distribution, correspond to the ones 
using the minimum time necessary to transit within a lobe touching a certain number of 



times the boundary to the other lobe, as illustrated in Fig. [9p. In Fig. [10} we illustrate the 
type of trajectories with major frequency in the distribution of Tperm- In this interesting 
case, we observe the presence of two kind of partial trajectories that go across the lobe 
entering through one boundary and going out through the same boundary; just one of them 
almost arrives to the second boundary. The other class of trajectories has a non negligible 
probability of being followed by an almost identical trajectory under reflection when the 



atom crosses to the other lobe. This is clearly illustrated in Fig. 10 d where consecutive 
partial trajectories of this kind lead to a partial trajectory very similar to the quasiperiodic 
case shown in Fig. |3b. 
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Figure 6: Power spectra of the coordinate x{t) for the three chaotic trajectories shown in 
Figure [5] The broad band spectra is a signature of chaotic motion. Only the low frequency 
end of the spectra is shown. The frequency is measured in cycles per unit time F"^. 
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V. CONCLUSIONS 



We have numerically shown that optical lattices with transverse structure allow the design 
of potentials yielding non trivial dynamics. For chaotic trajectories, like the one shown in 
Figure |5]3, the atom has a non zero probability to be found at any point within its accessible 
region of phase space. This behavior is similar to that expected for the quantum description 
of the atom motion. In this case the atom wave function will necessarily lead to a position 
probability density that mimics the optical structure, being non zero in any of the lobes. 
We expect that for ultra cold atoms chaos-assisted bidimensional tunneling will take place. 

Under the proposed set-up, one can notice that the fraction of trapped atoms with chaotic 
trajectories diminishes as the irradiance is lowered. Fig. [2j It is also important to emphasize 
that, for intermediate values of the irradiance, most of the trapped chaotic trajectories 
involve more than one lobe and that the distribution of time permanencies within a lobe 
is highly structured and with a long tail. Finally, for high irradiances, the proportion of 
chaotic trajectories saturates and there is a high presence of quasiperiodic trajectories. In 
this last regime, the atom is highly confined in the 2;-direction and tends to move inside a 
single optical potential lobe. Since the detuning is large, dissipative effects are expected to 
be small, thus, in the high intensity regime, the system dynamics is very similar to that of 
two-dimensional billiards with a topology determined by the parabolic symmetry of the light 
beam. It is known [26] that this kind of dynamical systems have a large set of quasiperiodic 
solutions. We consider this the key to understand the saturation effect illustrated in Fig. |2} 
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Figure 7: Trajectory in the velocity Vy versus position y phase space plane for the path 4c. 
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Notice that the use of a propagation invariant EM wave is not mandatory for the qual- 
itative results we found. The use of other beams like the Gaussian parabolic are expected 
to yield similar results for the transverse motion of atoms in the focusing plane. 

We also assumed anisotropic initial velocities. If the mean initial z-velocity were com- 
parable with the transverse velocities, one would expect that the atoms motion would take 
place in multiple quasi bidimensional configuration spaces since the atom would visit several 
of the sinusoidal potential wells in the z-direction. However, the projection of the trajec- 
tories in a given transversal plane would be expected to be qualitatively similar to those 
described above. 

For the geometries we are studying, blue detuned beams, due to the numerous routes 
of escape provided by the dark light zones, in general do not lead to trapped trajectories. 
Nevertheless this scattering process deserves a detailed study on its own, both classically 
and quantum mechanically. 




logio(Arr) 

Figure 8: Log-log plot illustrating the statistics of the time an atom expends in an optical 

lobe before jumping to another optical lobe for a chaotic trajectory. The sampling was 
taken for the trajectory illustrated in Fig|5)D after an analysis that considered a total time 
Tp = 2.5 X lO'^F"^ in which 82951 transitions took place. The frequency corresponds to the 

number of events on which the atom expended a time [T, T + 50F^^] within a lobe. The 
long tail is a signature of a statistics where the momenta is not well defined. In this figure, 

the logarithm of zero has been mapped to zero. 
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